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Abstract. In the absence of GMB precision measurements, a Taylor expansion has often been 
invoked to parametrize the Hubble flow function during inflation. The standard “horizon 
flow” procedure implicitly relies on this assumption. However, the recent Planck results 
indicate a strong preference for plateau inflation, which suggests the use of Pade approximants 
instead. We propose a novel method that provides analytic solutions of the flow equations 
for a given parametrization of the Hubble function. This method is illustrated in the Taylor 
and Pade cases, for low order expansions. We then present the results of a full numerical 
treatment scanning larger order expansions, and compare these parametrizations in terms 
of convergence, prior dependence, predictivity and compatibility with the data. Finally, we 
highlight the implications for potential reconstruction methods. 
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1 Introduction 

Inflation is one of the leading paradigms for explaining the physical conditions that prevailed 
in the very early Universe [1-6]. It consists in a phase of accelerated expansion that solves 
the standard hot Big Bang model problems, and provides a causal mechanism for generating 
inhomogeneities on cosmological scales [7-12]. The most recent Planck measurements [13-15] 
of the Cosmic Microwave Background (CMB) indicate that these cosmological perturbations 
are almost scale invariant, with undetected level of non-Gaussianities and isocurvature com¬ 
ponents. At this stage, the full set of observations can therefore be accounted for in the 
minimal setup, where inflation is driven by a single scalar inflaton field (j) with canonical 
kinetic term, minimally coupled to gravity and evolving in some potential V{4>)- 

However, since the inflationary mechanism is supposed to take place at very high ener¬ 
gies, in a regime where particle physics is not known and has not been tested in accelerators, 
the physical nature of the inflaton and its relation with the standard model of particle physics 
and its extensions remain elusive. The only condition on V is that it should be sufficiently 
flat to support inflation, but otherwise the multitude of inflaton candidates (with associated 
potentials) makes the theory as a whole hardly tractable, unless one restricts to a specific 
model or scan them all one by one [16-19]. 

Another strategy consists in developing model independent approaches and in studying 
generic parametrizations of the inflationary dynamics (see, for instance, Refs. [20-25]). Given 
the large volume in parameter space left by the first CMB anisotropy measurements [26-28] , 
Hoffman, Turner and Kinney [29, 30] proposed to study “generic” Hubble flow dynamics by 
performing a Monte Carlo analysis over a set of Hubble functions. This procedure, dubbed 
“horizon flow”, was shown to yield the “typical” predictions r ~ 0 or r ~ 16(1 — ns)/3, 
where ns is the scalar spectral index and r is the tensor-to-scalar ratio. When r ~ 0, it 
was found that generically, ns < 0.85 or ns > 1. However, Liddle pointed out [31] that this 
parametrization implicitly relies on Taylor expanding the Hubble function to some order, 
implying some specific potentials. In Refs. [32-34], a more detailed analysis of the problem 
revealed that the phenomenological class of inflationary potentials sampled by the procedure 
is indeed responsible for these predictions. In particular, it should not come as a surprise 
that potential reconstruction methods based on the horizon flow procedure naturally give 
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rise to chaotic-like potentials, despite the data. In some sense, those are already oversampled 
in the prior of the method. 

Phenomenological parametrizations of inflation always contain some bias towards a 
specific class of dynamics and it is therefore often tricky to blindly use them in order to 
obtain physical information from, say, observational data. This is why in this paper, we 
reverse the problem and ask the following question. Given what we observationally know, 
what is the most sensible parametrization of inflation that incorporates current observational 
constraints? What prescriptions should be used when studying other aspects of the early 
Universe, where one needs to model inflationary dynamics phenomenologically? 

The Planck satellite results have attracted attention to a particular set of inflationary 
scenarios [17, 35], plateau inflation. One possibility to include this piece of information in the 
way we parametrize inflation is to make use of Pade approximants of the Hubble function 
instead of Taylor expansions. The rest of this paper is therefore organized as follows. In 
section 2, we set the main notations and propose a new method to analytically integrate the 
flow equations. This relies on making use of the field redefinition invariance of the problem, 
in order to identify integration constants a la Noether. In section 3, we apply this method at 
low order to Taylor and Pade expansions of the Hubble function, and compare both results. 
(In appendix A, we also give the results for an inverse Taylor expansion, in order to further 
illustrate how our method works in practice.) Then, in section 4, we provide a numerical 
analysis of these two parametrizations, when generalized to higher orders. In particular, we 
compare the way they sample inflationary parameter space as a function of the truncation 
order, discuss their prior dependence and comment on their agreement (or lack thereof) with 
observations. Finally, in section 5, we recap our main results and draw a few conclusions, 
including comments on potential reconstruction with this method. 


2 Hubble Flow Dynamics 


The horizon flow formalism relies on the introduction of a set of “flow parameters” charac¬ 
terizing the way the Hubble scale evolves in time. There are several possible sets of such 
parameters. For example, let us consider the so-called “Hubble flow parameters”, defined by 
the flow equations [36, 37] 


din Cj 
^ diV 


( 2 . 1 ) 


The hierarchy is started at cq = Hin/H, where = In a is the number of e-folds and increases 
as time proceeds. Since ei = —H/H^ = 1 — d/{aH‘^), inflation (a > 0) takes place provided 
ei < 1. These parameters can be written in terms of the Hubble function H{(j)) and its 
derivatives. For the lowest orders, one has 


ei 


£3 



fH'V H"' H"' 

\h) 


) 


( 2 . 2 ) 


Throughout this paper we set the reduced Planck mass Mpi = 1, dots mean derivating with 
respect to cosmic time and primes refer to derivating with respect to the inflaton field (j). 

The horizon flow strategy rests on solving a truncated hierarchy of flow equations for 
a given set of flow parameters. Since these flow parameters can always be written in terms 
of the function and its derivatives, as in Eqs. (2.2), this procedure thus relies on 

the assumption that some combination of H and its derivatives, corresponding to the hrst 
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vanishing flow parameter, is zero. Interpreted as a differential equation for the H((p) function, 
this means that H{(j)) is parametrized in a certain manner, involving a finite number of 
constant free parameters. 

For example, if the Hubble flow hierarchy is truncated at some order M, i.e. if one 
assumes e; = 0 for Z > M, then em is constant, and H oc exp(ai exp(o 2 • • • exjp{aMN)' ")) 
where the exponential function is composed M times. As another example, if one makes use 
of the parameters, widely used in the horizon flow literature and defined as [30] 




^ W d(/)Hi ’ 


(2.3) 


truncating the hierarchy at order M means that = 0, hence H{(p) has a 

polynomial form of degree M. In general, one can see that truncating a specific flow hierarchy 
always boils down to parametrizing the Hubble function in a specific manner. 

Conversely, to any parametrization of the Hubble function, one can associate a specific 
dynamical system. Let H {cp, ai, a 2 , ■ ■ ■, an) be a given parametrization, where the coefficients 
Uj stem from some (e.g. Taylor or Fade ) expansion truncated at some order n. The n + 1 
first derivatives of this function with respect to the inflaton field (p can be calculated, and 
one can invert the system to extract the n + 2 variables {0, oi, 02 ,..., a^, | in terms of 

|H, H', H" ,..., }. Of particular interest is the last entry of this solution, which relates 

the (n + l)**^ derivative of the Hubble function to the lower order ones. In terms of the 

Hubble flow hierarchy, this means that €n+i can be expressed in terms of the n first Hubble 
flow parameters only. The flow equations (2.1), for 1 < i < n, thus form a closed dynamical 
system. It is important to stress that all physical input resides in this truncation: how En+i 
is expressed in terms of all preceding flow parameters fully determines the dynamical system 
and hence the inflationary predictions. 

Moreover, the flow dynamics is insensitive to the actual value of the inflaton field cp and 
hence the transformation cp —)■ (p-\-S(p leaves this system invariant. For the expansions that we 
consider, indeed, the functional form of H((p) does not change under this shift, implying that 
there is a degeneracy in the parameters. Amongst the n parameters, one combination can 
therefore be absorbed by the shift transformation, while the remaining n — 1 combinations are 
invariant. The later are therefore constants of motion, and the space of inflationary solutions 
has dimension n — 1. 

In the following section, we show how these constants of motion can be derived in 
practice, and we apply this method to second order Hubble and Fade expansions of the 
Hubble function. In these cases, a single constant of motion will be obtained. Therefore, at 
fixed number of e-folds AA^* between the Hubble exit time of the CMB pivot scale and the 
end of inflation, a one-to-one relation between ns and r will be obtained. 

A last remark is in order. Though the systems will be solved exactly and independently 
of the slow-roll approximation, in practice, the scalar spectral index and the tensor-to-scalar 
ratio will be calculated from the flow parameters at Hubble exit time thanks to the rela¬ 
tions [38, 39] 


Us — 1 — 2ei* — €2* — — (2C -|- 3)ei*e2* — Cg2*^3* i — 16ei* (1 -|- Ct2*) j (^-4) 

which are valid at second order in slow roll. Here, C = 7E-|-ln2 — 2~ —0.7296, 7 e being 
the Euler constant. In the Flanck preferred regions, the above expressions are valid and can 
therefore be used to compare predictions with observations. 
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3 Analytical Integration of the Hubble Flow 


In this section, we appiy the method sketched in section 2 to two toy cases: a Tayior expansion 
of the Hubbie function at quadratic order and a Fade expansion at iinear order. We obtain 
anaiyticai expressions for the inflationary trajectories in the parameter space (ei, € 2 ), as weii 
as for the number of e-foids reaiized aiong these trajectories. Finaiiy, we dispiay in each case 
the corresponding vaiues of ns and r. 


3.1 Taylor Expansion 

We first illustrate our method by considering the well-known case of a quadratic Hubble 
function: 

H = Ho{l + a^ + b^‘^) . (3.1) 

This case was already solved in Ref. [34] but employing a different approach. Making use of 
Eqs. (2.2), the two first flow parameters are given by 


/ a + Ibcj) \ ^ / a + 2b4> \ ^ 26 

V 1 -|- a 0 -|- bcf)^ J ’ V 1 -|- a(/> -|- bcfP' ) 1 -|- -|- bcfP' 


(3.2) 


Let us derive the constant of motion. Under the inflaton shift transformation (f) ^ (j) + 64>, 
the functional form of Eq. (3.1) remains unchanged if the coefficients of the expansion change 
according to 


(1 —y 


a + 2b6(l) 

1 -|- abcj) + bbcf)'^ ’ 


6 

l + a6^ + bdcj)"^ ’ 


(3.3) 


where Hq has also to be rescaled according to Hq —)• Hq (l -|- a54> -|- bbcj)'^). If one moves to the 
specific gauge where a vanishes, i.e. if one takes 64> = —a/{2b), then the shifted 6 coefficient, 
6 V (6 — a^/4), is gauge invariant. This implies that the following combination is invariant 
under the inflaton shift^: 


326^ 

^ “ a2 -46 


( 2 ei - 62 )^ 
£2 - £1 


(3.4) 


where the second equality has been obtained from Eq. (3.2). 

As pointed out in section 2 , a given parametrization of the Hubble function can be 
translated into a specific dynamical system. For the case at hand, since Eq. (3.1) implies 
that H"' = 0, Eqs. (2.2) give rise to 


£3 = 3ei 



(3.5) 


This truncates the infinite set of flow equations (2.1) for all the e* into a set of two differential 
equations for the first two flow parameters: 


d£i 

diV 


= £i £2 , 


de 2 

diV 


—26^ -|- 3eie2 . 


(3.6) 


This dynamical system generates a flow through the two dimensional space ( 61 , 62 ), which 
is displayed in the left panel of Fig. 1. Different trajectories can be labeled by different 
values of the invariant parameter 7 . It can easily be checked that Eqs. (3.6) indeed leave this 
particular combination invariant. 

^The case with a? — Ah or, equivalently, ei = 62, is singular and needs to be treated separately. It is 
straightforward to show that, in this case, one simply has ei* = £2* = 1/(1 + AW). 
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Figure 1. Second order Taylor expansion for the Hubble function. Left panel: flow lines of the 
system (3.6) in the plane (ei,e 2 ). The arrows indicate in which direction inflation proceeds. The blue 
dashed line corresponds to ei = 1 where inflation stops. The three regions I, II and III refer to the 
discussion in the main text. Right panel: Observational predictions in the (ns.,r) plane, compared 
with the Planck 2015 la and 2a contours. The green lines stand for the values o/ris and r computed 
40 e-folds before the end of inflation (dashed line), 50 e-folds (solid line) and 60 e-folds (dotted line). 
The grey segment at the bottom right stand for the fixed points (ei = 0 ,e 2 < 0). 


Moreover, by inverting Eq. (3.4), one can relate one of the remaining two slow-roll 
parameters to the other: 

62 = 2 ei + I + IV72 + 4761, ( 3 . 7 ) 

where ^ = ±1 = sign [( 2 ei — 62 ) 62(61 — 62 )] and the argument of the square root is always 
positive. Inserting Eq. (3.7) into the first of Eqs. (3.6) then leads to a first order differential 
equation for 61 (A^) that can be solved, and one obtains 


AN,, = N (cgend) - N ( 61 *) , 


iV(6i) 


2 

7 + ^ 1/72 + 4761 


+ -log 
7 


-^ 72 + 4761 - 7 
1/72 +4761 +7 


(3.8) 


where 6 i^end = 1 (we only consider cases where inflation has a graceful exit). This expression 
can be inverted to yield 61 as a function of N, 

-Wx ^/-Iif7(2e,-e2)<0 

‘ + * \o if 7 ( 2 ei -£ 2 ) > 0 ’ '' 

where x determines which branch of the Lambert function is to be used. Note that the 
sign of the combination 26 i — 62 appearing in the definition of x does not change during 
inflation, as the Hubble flow equations imply that d(26i — e 2 )/dN = 6i(26i — 62 ). Therefore, 
the branch of the Lambert function does not change during inflation. One can check that 
the above formula matches Eqs. (66-70) of Ref. [34] where it was first derived. 

Let us now discuss the structure of the phase space diagram plotted in the left panel of 
Fig. 1. According to the type of Hubble function one is dealing with, three possibilities must 
be distinguished: 
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• Firstly, if €2 > 2 ei, then ei vanishes in the far past while 
€2 takes a positive value. During inflation, both mono- 
tonically increase. This corresponds to a Hubble function ^ 
that has an inverted parabolic profile whose maximum is 
positive (region I). In this case, 7 > 0 and ^ = + 1 . 

4 > 

• Secondly, if ei < 62 < 2 ei, then both slow-roll parameters 
are vanishing in the far past and monotonically increasing 
during inflation. The corresponding Hubble function has a ^ 
parabolic profile with a negative minimum (region H). In 
this case, 7 > 0 and ^ = —1. 

4 > 

• Thirdly and lastly, if €2 < ei, both slow-roll parameters are 
again vanishing in the far past. However, during inflation, 
ei reaches a maximum, and then decreases back to zero, 
while 62 asymptotes to a negative value in the future. This ^ 
stems from a Hubble function with a positive minimum 
(region HI). In this case, 7 < 0 , while ^ = +1 before ei 
crosses its maximum and ^ = — 1 afterwards. 

These three regions are shaded with different colors in the left panel of Fig. 1 . One should 
note that thanks to the conservation of the sign of 7 defined in Eq. ( 3 . 4 ), a given inflationary 
trajectory never changes region. Amongst the third category, one can distinguish two cases. 
If the maximum value of ei is smaller than one, inflation never ends and reaches one of the 
fixed points (ei = 0,62 < 0 ). If; on the other hand, the maximum value of 61 is larger than 
one, and if one starts inflating with 62 > 0 , then inflation ends naturally when 61 = 1 . This 
happens when 7 < — 4 . 

Combining Eqs. (3.9), (3.7) and (2.4), the inflationary predictions of this class of models 
can be obtained and are displayed in the right panel of Fig. 1 for 40 < AiV* < 60. One should 
note that Eqs. (2.4) also make use of 63, but 63 is related to 61 and 62 thanks to Eq. (3.5). 
For “large-field” scenarios (region H), r is too large, and the model asymptotes the line 
r = 16(1 — res)/3 mentioned in section 1 and commented on in Ref. [34], which separates 
regions H and HI. For “hilltop” or “small-held” scenarios (region I), r is small, but ris is 
generically too red. When interpolating between these two cases, there is a small range of 
models for which r ~ 0.1 and the spectral index ns has marginally the right value. However, 
one can check that this corresponds to very hne-tuned initial values of the how parameters 
(or, equivalently, values of 7). Moreover, we will hnd in section 4 that higher-order terms 
change this result signihcantly. 





3.2 Fade Expansion 

Let us then further illustrate our method by considering the case of a hrst order Fade ex¬ 
pansion, 


H 

Ho 


1 + acf) 
l + hS 


(3.10) 


This case has not been considered in the literature before and provides a simple implemen¬ 
tation of the idea of “plateau inhation”. The hrst two slow-roll parameters can be obtained 
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from Eqs. (2.2), and one has 


2 (a — 6 )^ 

(1 + acj))'^ (1 + ’ 


4 (a — 6 ) (a + 6 + 2abcj)) 

(1 + a 4 >)‘^ (1 + 60 )^ 


(3.11) 


Here, we follow exactly the same approach as the one used for the Taylor expansion in 
section 3.1. For instance, under shift transformations 4> —?• (j) + 5(p, the functional form (3.10) 
is unchanged provided 


Cl —y 


a 

1 + a6(j) ’ 


b 

l + b6<j) ’ 


(3.12) 


where Hq is also rescaled according to Hq —)> Ho{l + a64 >)/(I + bScp). By moving to the gauge 
where the constant term in the numerator of Eq. (3.10) vanishes, i.e. 54) = —Ija, the b 
coefficient becomes 6/(1 — 6 /a), which is therefore gauge invariant. This implies that 


16\/2a6 62 — 4e^ 

|o- 6 | ^ ^3/2 


(3.13) 


is a constant of motion and can be used to label the different trajectories. 

Let ns recall that a given parametrization for the Hubble function can always be cast 
in a single dynamical system in the flow parameters space. For the present case, making use 
of the same procedure as before, Eq. (3.10) implies that H'" = /{2H'), and Eqs. (2.2) 

give rise to 

63 = — + -62 . (3-14) 

62 4 


Again, this truncates the dynamical system to a closed set of differential equations for (ci, 62 ), 
given by 


d 6 i 

dJV - '1'^ • 


d£2 _ 2,32 
dA^ “ ^ 4^2 ■ 


(3.15) 


In particular, one can check that the combination 7 defined in Eq. (3.13) is left invariant. The 
equation for 61 is the same as in Eq. (3.6), since it just defines 62 . However, the equation for 
62 is different. In general indeed, only the flow eqnation for the last flow parameter encodes 
the physical information abont the model, and propagates back to yield a specific dynamics 
for all flow parameters. The integrated flow lines of the above system are displayed in the 
left panel of Fig. 2. 

Let us now see how this system can be integrated analytically. By inverting Eq. (3.13), 
one can express 62 as a function of 61 , 


£2 



(3.16) 


where ^ = ±1 = sign( 62 ) changes when 61 crosses its minimum value. As before, inserting 
Eq. (3.16) into Eqs. (3.15) yields a first order differential equation for 61 (A^) that can be 
solved, and one obtains 


AAi* = A(6i,end) - A^(ei*), 


36^ 72 


(3.17) 


where again 6 i^end = 1- Contrary to the result obtained in section 3.1, this expression cannot 
be inverted analytically. 
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Figure 2. First order Fade expansion for the Hubble function. Left panel: flow lines of the sys¬ 
tem (3.15) in the plane (ei,e 2 ). The arrows indicate in which direction inflation proceeds. The blue 
dashed line corresponds to ei = 1 where inflation stops. The three regions I, II and III refer to the 
discussion in the main text. Right panel: Observational predictions in the (ns,r) plane, compared 
with the Planck 2015 la and 2a contours. The red lines stand for the values of and r computed 
40 e-folds before the end of inflation (dashed line), 50 e-folds (solid line) and 60 e-folds (dotted line). 
The grey dot at the bottom right stands for the fixed point (ei = 0,62 =0). 


Let us now discuss the structure of the phase space diagram ploted in the left panel of 
Fig. 2. According to the type of Hubble function, three possibilities must again be distin¬ 
guished: 


• Firstly, if €2 > 2ei, both ei and €2 increase as inflation 
proceeds, from the (repulsive) fixed point ei = 62 = 0 
reached in the infinite past. Inflation ends naturally when 
ei = I. This implements the idea of “plateau inflation” 
where the Hubble function is concave with a non vanishing 
plateau where inflation proceeds (region I). In this case, 
7 > 0 and ^ = + 1 . 

• Secondly, if —2ei < 62 < 2ei, ei = 1 is reached both in 
the past and in the future. In between, a finite period of 
inflation takes place where 62 increases and ei goes through 
a minimum. The corresponding Hubble function is convex 
and vanishes before the plateau is reached (region H). 
In this case, 7 < 0 while ^ = — 1 before ei reaches its 
minimum, and ^ = —I afterwards. 

• Thirdly and lastly, if €2 < —2ei, ei decreases as inflation 
proceeds, while €2 increases. The (attractive) fixed point 
ei = 62 = 0 is reached in the asymptotic future. This 
corresponds to a convex Hubble function for which the 
plateau is positive (region HI). In this case, 7 > 0 and 
^ = - 1 . 
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In particular, making use of Eq. (3.17), one can check that an inhnite number of e-folds can 
be realized in cases I and III. However, in case II, only a finite number of e-folds can be 
obtained. Parametrizing a given trajectory within region II by e 2 ,end = \/7 + the value of 
the second flow parameter at the end of inflation, this number is given by 


Nr, 


8 

2 f^2,end 


1 2 — 

^2,end 
{4,end - 4 ) 


(3.18) 


As expected, this number vanishes when 62 ,end approaches 0 and diverges when 62 ,end ap¬ 
proaches 2 . 

Combining Eqs. (3.17), (3.16) and (2.4), the inflationary predictions of this class of 
models can be obtained and are displayed in the right panel of Eig. 2 for 40 < < 60. 

Again, in Eqs. (2.4), £3 is related to ei and 62 thanks to Eq. (3.14). When inflation proceeds 
in region I, in the limit e 2 end I 5 one recovers the “typical” predictions of plateau inflation 
where r is small and Uq is in good agreement with the observational constraints. In this limit, 
Eq. (3.16) gives rise to £2 — , and one has 


/ 4 

UvrAivJ 


£ 2 * 




4 

3AW ■ 


(3.19) 


This translates into ns ~ 1 — 4/(3AA*) and r ~ AA* <C 1, which is what one would 
expect from a plateau inflation model with a !/(/> fall-off [23, 24]. 

Finally, let us note that this regime is interesting because £ 2 end 1 means that the last 
stage of the inflationary phase is realized far away from slow roll (let us recall that, here, the 
inflationary dynamics is solved without resorting to the slow-roll approximation). However, 
the number of e-folds realized between the time when £2 = 1 and the end if inflation when 
£1 = 1 can be calculated thanks to Eq. (3.17), and in the limit where £ 2 end 1, one obtains 
4/3 e-folds. This is why, A A* e-folds before the end of inflation, slow roll is well valid and 
the system gives rise to predictions that are in good agreement with observations. 


3.3 Comparison 

In order to summarize the analysis of the two toy models discussed in the present section, in 
Fig. 3, we have superimposed their predictions in the (ns,r) plane, for a few fixed values of 
A A* G [40, 60]. One can see that Taylor and Fade lines are tangential, along the £3 = £2 = 2£i 
line, which is associated to the model H/Hq = cj). This should not come as a surprise for the 
following reason. After a suitable gauge transformation, Eq. (3.1) can be cast in the form 


H 

Ho 



(3.20) 


where 7 has been dehned in Eq. (3.4). Similarly, Eq. (3.10) can be cast in the form 


H _ (j) 


(3.21) 


where 7 has been defined in Eq. (3.13). As a consequence, a linear Hubble is a special case 
of both parametrizations, corresponding to 7 = 0. However, obviously, the way 7 modifies 
this linear function is different for both parametrizations. 
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Figure 3. Compared predictions of the Taylor model (3.1) (green) and the Fade model (3.10) (red), 
for 40 < AiV* < 60, in the (ns,r) plane. The black lines are the la and 2a contours of Planck 2015. 
The grey solid line stands for the model H/Hq = <)>, which is a special case of both parametrizations 
and for which £3 = €2 = 2ei. 


4 Numerical Integration of the Hnbble Flow 

The above results indicate that inflationary dynamics is better parametrized by Fade expan¬ 
sions of the Hubble function rather than Taylor expansions. However, one might worry that 
this statement relies on the low truncation order we have worked with. This is why in this 
section, we generalize our approach by including higher order terms numerically. This also 
allows us to investigate a crucial aspect of flow parametrizations, namely the dependence (or 
lack thereof) of the results on the prior choice for the parameters of the expansion. Notice 
also that the analytical method developed in section 3 can in principle be used to deal with 
arbitrarily large order expansions, however, one would not gain much by displaying the cor¬ 
responding cumbersome formulas. This is why, here, we directly compute the predictions of 
the models we study, which consist in Hubble functions of the form 

H{4>) = H [m, n] {ct>) = • (4.1) 

Ez=o di(t^ 

In practice, we consider orders [M, 0], which correspond to Taylor expansions of the Hubble 
function, and orders [M, M], which correspond to Hubble functions that asymptotes to a 
non-vanishing plateau at large-held values. We study different values of M in order to test 
the robustness of the predictions under increasing the order of truncation and we report the 
results below. 

When the Hubble function is described by a Taylor series, we start our exploration 
at (/>ex = 0. This is because, in principle, the Taylor expansion is dehned around (/> = 0, 
and might not always converge far from this value. A Fade approximant is a simultaneous 
expansion around 0 = 0 and 0 = 00 , which allows us to either start the exploration of the 
Fade Hubble function at 0ex = 0 or 0ex = 00 . This last case corresponds to plateau inhation 
and we expect it to be in best agreement with the data. A hrst how algorithm was designed 
in Ref. [30], and here we use the modihcation proposed in Ref. [33] that allows us to include 
generic Hubble functions. It proceeds according to the following steps: 
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(i) Draw the parameters Cj~ and appearing in Eq. (4.1) according to some prior distri¬ 
bution (see below). 

(ii) The Hubble function being specified, calculate ei{4>) through Eq. (2.2). 

(hi) If ei(0ex) > 1, restart from (i). Otherwise go to (iv). 

(iv) Calculate the set of values {(/)o} where ei = 0 and where ei = 1. 

(v) If > 0, inflation proceeds at decreasing values of (p. Amongst the elements 

of {(/>o} and {4>i}, find which value is the closest to (/>ex while being smaller than (^ex¬ 
it this value belongs to {(/>o}, a fixed point is reached and we chose not to consider 
such trajectories in the present analysis since the results then depend on the value of (p 
at which inflation terminates, which has to be determined by other physical processes 
than slow-roll violation. If, on the other hand, it belongs to identify ())end with 

this value. Apply a similar procedure if ff'(^ex) < 0. 

(vi) Calculate the value of (p, A A* = 50 e-folds prior to ei = 1, by integrating 

/■'/'end A A 

AN, = sign(H') / ^ (4.2) 

v4ei 

and inverting the result. If a sufficient number of e-folds cannot be obtained, go back 
to (i). 

(vii) Calculate the slow-roll parameters ei*, € 2 * and es* making use of Eqs. ( 2 . 2 ), and the 
values for ns and r thanks to Eqs. (2.4). 

In practice, we iterate this procedure until we obtain 10® successful realizations. Note that 
the results presented below have also been derived with different values of AN, , and that we 
have checked that our conclusions remain unchanged. 

In step (i), the coefficients of the expansion are drawn according to some priors that we 
now specify. We studied two classes of priors. The first one consists in drawing all coefficients 
Cfc and dk from a flat distribution between [—p/q^,p/q^], where/? and q are two fixed numbers. 
In the following, it is referred to as the “power-law” priors. In the second class of priors, Ck 
and dk are drawn from flat distributions in the interval [—pfkjPfk]^ where p is a constant 
and the set {fk} is defined such that if Ck = {—l)^dk = fk for all k, all M derivatives of the 
Hubble function ai (p = Q are 1. This gives rise to 

Taylor: fk = -j^, Fade: /^ =( 4 . 3 ) 

In what follows, this prescription is referred to as the “binomial priors”. 

The original procedure proposed by Kinney [30] relies on drawing Ck G {\/2/{{2q)^{k + 
!)!))[—p,p]/[ 0 , with do = 1 and dk = 0 ior k > 1 , where [—p,p]/[ 0 , stands 

for the ratio of two numbers, one drawn in the interval [—p,p] and the other one drawn in 
[0,1] and taken to the {k — l)/2 power. This corresponds to taking power-law priors on the 
values of the A flow parameters, defined in Eq. (2.3), at cpex = 0. The la contour obtained 
with this prescription is shown in the left panel of Fig. 4. It is widely spread, and since the 
2a contour is even more spread and difficult to resolve statistically, we do not display it. 
From the Planck contours, one can see that the vast majority of trajectories generated by 
this procedure are observationally excluded. This can be further quantified by calculating 
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Figure 4. Left panel: one sigma contour of the predictions generated by a Taylor expansion of the 
Hubble function, with priors on the coefficients of the expansion matching the procedure originally 
proposed in Ref. [30]. Far from rig = 1 and r = 0, the slow-roll formulas (2.4) are not valid and 
we display the contour for indicative purpose only. The Planck 1 and 2a contours are shown for 
comparison. Right panel: prediction points for a Taylor expansion, truncated at different orders M, 
when the coefficients are drawn from the binomial prior distributions. 

the percentage of the points lying outside the Planck 2a contour. One obtains 0.2%, a small 
fraction indeed. We should stress that these results correspond to the standard “horizon 
flow” procedure as commonly used in the literature. They again motivate our search for 
alternative parametrizations. 

Unfortunately, other natural prior choices for the Taylor expansion coefficients, such as 
the power-law or the binomial priors introduced above, fail to converge as M increases. In 
other words, the obtained results have a large dependence on the order of truncation, an 
undesirable property. More precisely, in the power-law case, since the radius of convergence 
of the Taylor series is (f = q, the model converges only when q is sufficiently large. For such 
priors, the higher order terms are strongly suppressed and, in practice, the simple results 
of the second order truncation obtained in section 3.1 are recovered, and the higher order 
information is quickly lost. In the binomial case, the results corresponding to different orders 
of truncation are displayed in the right panel of Fig. 4, where it is obvious that the model 
does not converge. Therefore, it seems difficult to design alternative parametrizations relying 
on a Taylor expansion. 

Let us now consider Fade approximants. First of all, we have checked that with the two 
classes of priors proposed above, the results always converge (indeed, it is well known that 
Fade approximants have better convergence properties than Taylor expansions). In practice, 
we hnd that M = 6 is enough to reproduce all higher order results with a very good accuracy. 
The next question is how much the results depend on the class of priors. In Fig. 5, we have 
displayed the 1 and 2a contours obtained with the power-law and binomial priors, in the case 
where 4>ex = 0 (left panel) and cpex = oo (right panel). 

When (/>ex = 0, the binomial prior gives rise to wide spread contours (the 2a contour 
entirely lies outside the plot frame). They are consistent with the power-law results, but scan 
larger sets of inflationary trajectories. For this reason, the percentage of points inside the 
Flanck 2cr contour is smaller, 5% for the binomial distribution and 38% with the power-law 
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Figure 5. 1 and 2a contours for the two prior choices (power law and binomial), for a Fade expansion 
of the Hubble function. The left panel corresponds to = 0 and the right panel to (/)ex = oo. 

prior. 

When (pex = OO; we find that both priors give rise to rather narrow contours, in agree¬ 
ment with each other at the Icr level. In this case, large fractions of points lie inside the 
Planck 2a contour: 18% for the binomial prior choice and 90% when the power-law prior is 
used. These models correspond indeed to one’s intuitive representation of “plateau inflation”. 
Interestingly also, a lower bound on r is found, which means that this class of inflationary 
dynamics could in principal be ruled out by future experiments. 

This analysis thus reveals that up to a moderate prior dependence. Fade expansions 
of the Hubble function give rise to predictions in agreement with observations, and possess 
good convergence properties. Given what we observationally know, they seem better suited 
to parametrize inflationary dynamics than Taylor expansions, on which the standard horizon 
flow procedure rests. 

5 Conclusion 

Let us now summarize our main findings. Whenever inflation is parametrized by a truncated 
dynamical system for the flow parameters, it can equivalently be described by an expansion 
scheme for the Hubble function at some finite order. Conversely, any functional shape 

for the Hubble function (such as a Taylor expansion, a Fade expansion, or any other expansion 
involving a finite set of free coefficients) can be related to a single dynamical system in the 
flow parameters space. 

Making use of the shift symmetry (p —)> p + 64> of the problem, we have explained 
how constants of motion can be derived for such systems, and how their dynamics can be 
integrated. For illustrative purpose, we have applied this new method to the case of a second 
order Taylor expansion (section 3.1), a first order Fade expansion (section 3.2), and a second 
order inverse Taylor expansion (appendix A). For the second order Taylor case, we have found 
that generically, either r is too large and the famous horizon flow relation r ~ 16(1 — res)/3 
is recovered, or r is small enough but Uq is too red. For the first order Fade expansion on 
the other hand, the typical predictions ns ~ 1 — 4/(3AW) and r ~ AW <C 1 have been 
obtained, in good agreement with observations. 
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Figure 6. Left frame: la and 2a contours for the two Fade parametrization around </> = 0 and 
(j) = oo, both with power-law priors. The predictions of the first order Fade expansion, as obtained 
in section 3.2, are displayed too. Right frame: potential reconstruction for the the most likely 100 
trajectories (using the Flanck likelihood [If]), for the standard horizon flow procedure based on a Taylor 
expansion of the Hubble function, a Fade expansion around (f = 0, and a Fade expansion around 
(j) = oo. Darker colors stand for smaller likelihood of the predicted values for Us and r. The potentials 
are normalized so that the amplitude of the curvature power spectrum, Vc^ = 14/(247r^ei*Mpj) ~ 
2.203 X 10“®, is correctly obtained. 

We have then extended these results to higher order expansions by numerical means in 
section 4, and studied the dependence of the predictions on the priors chosen for the coef¬ 
ficients of the expansions. We have confirmed that Fade expansions of the Hubble function 
are more suited to parametrize inflation than Taylor expansions, since they show good con¬ 
vergence properties, mild prior dependence, and, most notably, much better agreement with 
observations. 

When using Fade approximants, we have distinguished the case where inflation pro¬ 
ceeds close to (/) = 0 and close to f = oo. These two prescriptions give rise to results 
that are compared in the left panel of Fig. 6, where we also display the first order Fade 
result of section 3.2. It is clear that, given observational constraints, they provide a better 
parametrization of inflationary dynamics than the usual horizon flow procedure. 

These results illustrate why “model-independent” parametrizations of inflation, such 
as expansion schemes for the Hubble function or truncated flow dynamical systems, always 
make non trivial assumptions about its dynamics. This is why, a priori, parametrizations 
yielding predictions that are in contradiction with the data cannot be used to infer physical 
information from them. 

As an example, in the right panel of Fig. 6, we show the best 100 potentials, obtained 
from the formula V = , for the three cases: Taylor expansion with the 

priors corresponding to the usual “horizon flow” procedure. Fade expansion around cj) = D, 
and Fade expansion around (/> = oo. From here, it is clear that different parametrizations 
sample inflation along different classes of potentials, even when restricted to the best possible 
realizations. Taylor expansion and Fade expansions around 4> = 0 seem to prefer inflationary 
potentials with a flat inflection point. Because this is a rather fine-tuned conhguration from a 
generic Taylor expansion, we understand why the standard horizon flow procedure, which is 
based on a generic Taylor expansion, produces trajectories that are most of the time excluded 
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by observations. On the contrary, Fade expansions around (j) = oo mostly samples plateau 
inflation, as expected. 

As a consequence, it is clear that “reconstructing” the potential with either of these 
parametrizations biases the result towards the class of potentials that it relies on. This 
is in essence a “prior effect”. However, phenomenological descriptions are still very useful 
to address a number of other issues in the Early Universe, where the background is effec¬ 
tively “sourcing” some physical effects. Therefore, the question becomes: “How can we best 
obtain and parametrize a class of inflationary trajectories that are in agreement with cur¬ 
rent observational constraints?” From a Bayesian perspective, the priors for analyzing the 
r^th generation of data come from the information provided by the n — 1**^ survey. In this 
respect, we have shown that after Planck, Fade expansions (or other types of expansion 
schemes implementing the plateau structure) should be preferred over the standard Taylor 
parametrizations of the Hubble flow dynamics. 
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A Inverse Taylor Expansion 


In order to further illustrate the method depicted in section 2, in this appendix, we apply it 
to the case where the inverse Hubble function is Taylor expanded at second order. 


H _ 1 

Hq 1 + a4) + 


(A.l) 


The two first slow-roll parameters can be read off from Eqs. (2.2), and one has 


2 (a -|- 2b(j))^ ^2b‘^(fP' -|- 2ab(f) + — 2b 

JTTa^+Wf ’ ~ (1 + a(j) + b(j)‘^f 


(A.2) 


Then, under shift transformations 4> ^ (f) + 64>, the functional form (A.l) is unchanged 
provided 


a + 2b6(l) b 

1 -|- a6(j) -|- ’ 1 + a6(j) + b5cjP‘ ’ 

where Hq is also rescaled according to Hq —)■ Hq /{1 + a 64 > + bScfP'). These gauge transfor¬ 
mations are, for obvious reasons, the same as for the second order Taylor case studied in 
section 3.1, which implies that 


326^ _ (261 + e2f 
46 -a? 61 -b 62 


(A.4) 


is a constant of motion and can be used to label the different trajectories. 

Let us now derive the dynamical system associated to this parametrization of the Hubble 
function. Making use of the same procedure as in section 3, we find that Eq. (A.l) implies 
that H'" = QHH'H"/H^, and Eqs. (2.2) then gives rise to 


£3 = -ei 


3 + 2- 
£2 


(A.5) 
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This truncates the dynamical system to a closed set of differential equations for (ei, 62 ), given 
by 


dei 

dN - ’ 


de 2 

div 


—ei (3e2 + 2ei) . 


(A.6) 


In particular, one can check that the combination 7 defined in Eq. (A. 4) is left invariant. 
The integrated flow lines of the above system are displayed in Fig. 7. 

Let us now see how this system can be integrated analytically. By inverting Eq. (A. 4), 
one can express 62 as a function of ei, 


62 


= -2ei + I + |\/7^-47ei 


(A.7) 


where^ ^ = ±1 = sign 


£ 2 ( 2 ei+e 2 ) 

ei+C 2 


As before, inserting Eq. (A.7) into Eqs. (A. 6 ) yields 


a first order differential equation for ei{N) that can be solved, and one obtains = 

A^(ei,end) - A^(ei*), where ei,end = 1 and 


iV(ei) = 


+ T^log 


— dyei — 7 27 


y/72 _ ry 


\/'f‘ - 4761 _ 7 


(A. 8 ) 


As in the case of the Taylor parametrization, Eq. (A. 8 ) can be inverted. 


61 (A) 


-4W^ (-e-7AAf-i) , _ [0 if 7(261 + 62 ) > 0 

[1 + (-e-7AiV-l)]2 ’ if^(26i+62) <0 


(A.9) 


where x determines the branch of the Lambert function W^. It is easy to show that x does 
not change along a given trajectory. 

Let us now discuss the structure of the phase space diagram plotted in Fig. 7. Accord¬ 
ing to the type of Hubble function one is dealing with, two possibilities must be distinguished: 


• Firstly, if 62 < — 61 , 61 decreases and 62 increases as 
inflation proceeds, reaching one of the (attractive) fixed 
point (61 = 0 ,62 < 0) in the asymptotic future. The 
corresponding Hubble function has a convex shape with a 
positive minimum (region I). In this case, 7 < 0. 

• Secondly, if 62 > — 61 , 62 decreases as inflation pro¬ 
ceeds while 61 first increases, crosses a maximum and 
then decreases, reaching the (attractive) fixed point 
(61 = 0, 62 = 0) in the asymptotic future. The correspond¬ 
ing Hubble function is concave and vanishes at infinity 
(region H). In this case, 7 > 0. 




One should note that thanks to the conservation of the sign of 7 defined in Eq. (A.4), a 
given inflationary trajectory never changes regions. Amongst the second category, one can 
distinguish two cases. If the maximum value of 61 is smaller than one, inflation never ends 
and reaches the fixed point (61 = 0, 62 = 0). If, on the other hand, the maximum value of 
61 is larger than one, and if one starts inflating with 62 > 0 , then inflation ends naturally 


^This sign does not change for trajectories for which ei always decreases with time (regions I in the 
classification introduced below). However, for trajectories for which ei first increases, reaches a maximum and 
then decreases, it is positive before reaching the maximum and negative afterwards (region II). 
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Figure 7. Inverse second order Taylor expansion for the Hubble function: flow lines of the sys¬ 
tem (A.6) in the plane (£ 1 , 62 ). The arrows indicate in which direction inflation proceeds. The blue 
dashed line corresponds to ei = 1 where inflation stops. The two regions I and II refer to the discussion 
in the main text. 


when ei = 1. This happens when 7 > 4. In this case, one can check that the function N (ei) 
defined in Eq. (A. 8 ) goes to infinity when ei goes to 0 which means that a sufficient number 
of e-folds can always be realized. 

However, this requires ei* to be sufficiently small. If one parametrizes a given trajectory 
within region II by £ 2 ,end = —2 -|- (7 + — 47 )/2, the value of the second flow parameter 

at the end of inflation, in the ei* <C 1 limit one has 


/ 2 -|- e2end 
V 1 + e2end 


2 

exp 


1 — (2 + e2end) AAI* 

1 + £2end 


(A.IO) 


Since 62 * > e 2 end) ^ 2 * < 1 implies that e 2 end < 1 hence ei* < 10 if one lets AAi* = 60. So 
essentially, r ~ 0 in these models. Making use of Eq. (A.7), one then has 


(2 + e2end)^ 
1 + e2end 


(A.ll) 


Because e 2 end > 0 in those branches where inflation ends naturally, this means that £ 2 * > 4, 
which is of course completely excluded by CMB observations. As a conclusion, the only 
trajectories compatible with observations are those that reach the fixed points (ci = 0 , £2 < 0 ) 
such that £2 gives the correct value of n^. However, these are rather fine-tuned situations 
that moreover require to invoke an extra mechanism to end inflation. 
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